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Abstract. The energy of the thermal flare plasma and the kinetic energy of the non-thermal electrons in 14 hard X-ray peaks 
from 9 medium-sized solar flares have been determined from RHESSI observations. The emissions have been carefully sep¬ 
arated in the spectrum. The turnover or cutoff in the low-energy distribution of electrons has been studied by simulation and 
fitting, yielding a reliable lower limit to the non-thermal energy. It remains the largest contribution to the error budget. Other 
effects, such as albedo, non-uniform target ionization, hot target, and cross-sections on the spectrum have been studied. The 
errors of the thermal energy are about equally as large. They are due to the estimate of the flare volume, the assumption of the 
filling factor, and energy losses. Within a flare, the non-thermal/thermal ratio increases with accumulation time, as expected 
from loss of thermal energy due to radiative cooling or heat conduction. Our analysis suggests that the thermal and non-thermal 
energies are of the same magnitude. This surprising result may be interpreted by an efficient conversion of non-thermal energy 
to hot flare plasma. 
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1. Introduction 

In the standard model of solar flares, a major part of the en¬ 
ergy is first released into energetic, non-thermal electrons and 
possibly ions. These particles, guided by magnetic held lines, 
may be lost in interplanetary space, but also precipitate into the 
lower corona or upper chromosphere where they loose their 
energy by Coulomb collisions with the denser medium. This 
energy is believed to heat up the ambient plasma to tens of mil¬ 
lions of degrees, which will rise and All the coronal loop. An 
interesting question is thus the relation between non-thermal 
and thermal energies. The difference between them may indi¬ 
cate conversion losses and/or other forms of primary energy 
release. 

How does total kinetic energy of electrons that precipi¬ 
tated compare to the thermal energy of the plasma? In the pre- 
RHESSI era and in RHESSI first results papers, this issue had 
been addressed (see e.g. de Jager et al. Saint-Hilaire and 
Benz 2002), with the result that the kinetic energy was often 
reported to be up to an order of magnitude or more than the 
observable thermal energy. Others (Gan et al. 2001) have man¬ 
aged to conclude just the opposite. 

The hot thermal plasma emits soft X-rays. They differ in 
the spectrum from the X-rays emitted by the non-thermal elec¬ 
tron bremsstrahlung at higher energies. Soft and hard X-rays 
can thus be used to determine the thermal and non-thermal en¬ 
ergies, respectively. In practice, however, many uncertainties 
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limit the precision of the energy determination. Unfortunately, 
the two emissions are usually cospatial and overlap in photon 
energies in the range from about 10 to 25 keV. This is further¬ 
more the range of the emission of the non-thermal photons that 
carry the information on the dominant lowest-energy part of 
non-thermal energies. The distinction between the two forms of 
energy requires a high spectral resolution in the critical range. 

An accurate derivation of non-thermal energy from the 
spectrum is not trivial. In previous work the observed hard X- 
ray spectrum was converted into electron energies assuming 
that they impinge onto a thick target. The energies were inte¬ 
grated starting at some assumed lower limit or at the crossover 
between thermal and non-thermal emissions in the observed 
photon spectrum. As the derived electron energy distribution is 
a power law with large index, the result depends much on the 
lower bound of integration. In addition, the photon spectrum is 
influenced by several effects. The non-thermal bremsstrahlung 
of a coronal source is reflected by the dense layers below, and 
thus appears to be brighter (Bai & Ramaty 1978; Alexander 
& Brown 2002). Precipitating electrons above a certain energy 
penetrate into the chromosphere where they lose energy by col¬ 
lisions with neutrals and are more efficient in bremsstrahlung 
radiation than in the completely ionized corona (Brown, 1973; 
Kontar et al. 2002). Additionally, the Sun is not a simple cold, 
thick target (Emslie 2003). Einally, the various approximations 
for electron cross-section used in the literature yield different 
values in particular at relativistic particle energies. 

On the other hand, the accurate determination of the ther¬ 
mal energy also poses problems. Thermal energies are best es- 
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timated from the thermal bremsstrahlung spectrum, yielding 
the plasma temperature and the emission measure. The volume 
of the thermal plasma must be estimated from the size of the 
source and an assumption on the hlling factor. The ambient 
plasma cools down by either heat conduction or radiative cool¬ 
ing (i.e. thermal bremsstrahlung of an optically thin plasma) 
(Porter and Klimchuck 1995; Aschwanden et al. 2001). Thus, 
thermal energy is lost as it is being measured. A reliable deter¬ 
mination of the errors is therefore as important as the hnal ratio 
between the two forms of energy. 

In this paper, we determine energy budgets for several 
flares observed by the Reuven Ramaty High Energy Solar 
Spectroscopic Imager (RHESSI) (Lin et al. 2002). The imag¬ 
ing capabilities, high spectral resolution and broad photon en¬ 
ergy coverage of RHESSI make this instrument ideal to de¬ 
termine the two energies with much higher precision and study 
their relation. Medium-sized flares (upper C and lower M class) 
have been selected to avoid photon pile-up, complicated source 
structure and other conundrums of larger flares. With a careful 
RHESSI analysis, some -but not all- of the uncertainties or as¬ 
sumptions used by previous authors may be removed, and more 
accurate results obtained. 

This paper will start with a section discussing some al¬ 
ready established facts of solar flare bremsstrahlung emissions 
and the various high-energy electron cross-sections available in 
the literature. The next section will deal with the method pro¬ 
posed to determine flare energies. Observational results from 
RHESSI will then be presented, and the ratios of cumulative 
non-thermal energy over thermal energy increase (during the 
same time interval) will be examined. 


2. Basic theory 

2.1. Thick-target bremsstrahlung emission 

An initial (injection) electron beam with distribution Fo{Eo) 
electrons s * keV ' passing through a dense -although still op¬ 
tically thin- plasma emits HXR radiation according to the fol¬ 
lowing formula, for a thick target (Brown 1971): 








Fig. 1. Synthetic photon spectra generated from perfect injec¬ 
tion electron power-law spectra of varying spectral index 6. 
Top row: The solid line were computed using the same non- 
relativistic cross-sections as Brown (1971). The dashed line 
were computed with the Haug (1997) cross-section and the 
Bethe-Bloch formula for energy loss. The bottom row displays 
the ratio between the two. 


where Ag is in electrons s“* keV"', is 1.44 for typical 
coronal abundances, kbh - ^ar^^mgC^ = 7.9 x 10“^^ cm^ 
keV, K - 2ne‘‘K - 2.6 x 10“'^ cm^ keV'^ for a fully ion¬ 
ized plasma, and B is the beta function [B{a,b) - 
(Tandberg-Hanssen and Emslie 1988). Numerically, 

hhickie) = 1.51 X 10-34 ^^*^ 

(d- l)(d-2) 

Figure ^ shows the difference when using more accurate 
relativistic cross-sections: The Haug (Haug 1997) differential 
bremsstrahlung cross-section and the full quantum relativistic 
Bethe-Bloch formula (Longair 1992) for energy loss: 

"(S) + 


^thick ~ 


AttD^ 


r r- .r-. eB(E,e) 


( 1 ) 


where Ithicki^) is the observed photon spectrum seen at distance 
D (assumed to be 1 AU in the following) from the site of emis¬ 
sion, in photons s ' cm“^ keV '. Isotropic emission is assumed 
throughout this work. Qb{E, e) is the bremsstrahlung differen¬ 
tial cross-section, for an electron of energy E emitting a pho¬ 
ton of energy e. QdE) is the energy-loss cross-section due to 
Coulomb collisions with the ambient plasma. 

Using non-relativistic cross-sections. Brown (1971) has 
demonstrated that there exist an analytical relationship between 
Ithick(e) and Fo(Eo), if Fo(Eo) is a power-law. If Fo(Eo) - 
AeE~^, then I,hick{e) - A^ ■ with: 

^ A. ^kbh B(d-2,1/2) 

AttD^ K (6-l)(6-2) 

■y — 6 - 1 


( 2 ) 

(3) 


where A is the usual non-relativistic Coulomb logarithm, y the 
Lorentz factor, and (3 — In its non-relativistic (NR) limit, 
the Bethe-Bloch formula has the ~ ^ dependence, up to elec¬ 
tron kinetic energies nearing the electron rest mass energy, af¬ 
ter which the dependence is ~ -g (Fig. |5}. We note here that 
the Bethe-Bloch cross-section is very close to what is used in 
RHESSI software. 

Other effects, such as photon back-scattering on the photo¬ 
sphere (Bai & Ramaty 1978; Alexander & Brown 2002), non- 
uniform target ionization (Brown 1973; Kontar et al. 2002), or 
a hot target (Emslie 2003) need to be considered. 

The non-thermal kinetic power of the injected electron 
beam may be computed: 

X OO 

E-Fo(E) dE (6) 

The introduction of some kind of a cutoff at low energies is 
necessary by the fact that the integral diverges at zero energy. 
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Fig. 2. Energy-loss cross-section, for the usual non-relativistic 
(NR) case and for the Bethe-Bloch formula. Notice also the 
slight dependence of the spectral index at relativistic energies 
with the ionization level. The spectral index at high energies 
is ~ 0.93 for a fully-ionized plasma, and ~ 0.87 for a neutral 
medium. 


Fig. 3. Top: Photon spectra and their spectral indices {bottom), 
computed from different injection electron power-laws of spec¬ 
tral indices 6, using the Brown (1971) cross-sections. The elec¬ 
tron spectra had all a 20 keV cutoff (solid line) or turnover 
dashed line energy. 


Assuming that Fq{Eq) remains a power-law below the thermal 
energy in the acceleration region is not physical. In the past, a 
sharp low-energy cutoff in Fo{Eo): 


FoiEo) = 


AeE^^ for Eo > Eco 
0 for Eq < Eco 


(7) 


was often assumed. This situation seems physically not realis¬ 
tic as such a conhguration leads to plasma instabilities. Such 
instabilities have growth rates typically of the order of the lo¬ 
cal plasma frequency, i.e. orders of magnitude shorter than the 
propagation time of the beam within the acceleration region. 
The (flat) turnover model seems physically closer to reality: 


Fo(£o) = 


AcEq^ forFo > E,o 
AeE;^ for Eq < E,o 


( 8 ) 


We note, however, that recent simulations of stochastic elec¬ 
tron acceleration, such as by Petrosian & Liu (2004) predict 
an energy distribution in the acceleration region still increas¬ 
ing below the turnover energy. Coulomb losses tend to lin¬ 
earize the low-energy part of an electron energy distribution. 
From their resulting photon spectra, such a linearization below 
the turnover would be extremely difficult to distinguish from 
a flat turnover, as it effects the least energetic photons, where 
the thermal emission usually dominates (see also the next para¬ 
graph). Thus, the flat turnover model may not be the true distri¬ 
bution, but yields rather an upper limit to the total non-thermal 
energy estimate. 

The photon spectra of the different model distributions of 
electrons in a beam impinging on a thick target have been 
calculated using the Brown (1971) model (i.e. using the non- 
relativistic Bethe-Heitler bremsstrahlung cross-section, and 
non-relativistic collisional losses). Both turnover and cutoff 
model lead to a photon spectrum that is rounded off at low en¬ 
ergies (see Fig. |3l. The spectral index of the turnover model 
is slightly larger than for the cutoff model. Note that it is not 


a constant value of ~ 1.5 as is sometimes assumed. Both tend 
asymptotically towards ~ 1.15 for all 6 as the photon energy ap¬ 
proaches zero. The usually observed superposition of a thermal 
component (or a full differential emission measure distribution) 
to the power-law spectrum further makes an observational dis¬ 
tinction in the spectrum exceedingly difficult. 

If Fq{Eq) has a cutoff with the shape dehned in Eq. Q. the 
non-thermal kinetic power contained in the beam of electrons 
is given by: 

Ecu,off - 2 . ( 9 ) 

On the other hand, if Fq{Eq) has the form of a turnover as de¬ 
scribed in Eq. (ID, 

P,urno.er^-^E;f^^{l + ^^'^ . ( 10 ) 

Figure|3]is an enlargement of the turnover region in the pho¬ 
ton spectrum. Note that the apparent photon turnover energy in 
Fig. a is below the electron cutoff/turnover energy (20 keV). 
This has been further investigated in Figs. |D and |D Synthetic 
photon spectra were produced from electron distributions with 
different power-law indices 6, using the Brown (1971) approxi¬ 
mation for easy comparison with analytical results. The spectra 
were then htted with a double power-law of spectral index 1.5 
below the break. The intersection of the two power-laws dehnes 
a photon turnover energy which is far below the electron cut¬ 
off energy Eco or turnover energy E,o. Fig. [^displays the ratio 
of the photon turnover energy e,o to either Eco or Eto- These 
ratios are smaller for the cutoff model than for the turnover 
model. If the photon turnover energy eto is used instead of Eco 
or Eto, the derived non-thermal power in the electron beam may 
be overestimated by more than an order of magnitude. 
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Fig. 4. Photon spectrum observed at 1 AU, produced by an elec¬ 
tron power-law distribution of spectral index 6=4, and a 50 keV 
electron flux of A 50 = 1.295 x 10^^ electrons keV^^. Solid 
line: No cutoff. Dotted line: Cutoff at 20 keV. Dashed line: 
Turnover at 20 keV. 



Electron spectral index 6 


Fig. 6. The ratio of the fitted turnover energy in the photon 
spectrum, to the cutoff resp. turnover energy in the electron 
beam distribution, as a function of the electron spectral index 
right of Ecojo- 



Fig. 5. Fitting a photon spectrum with power-laws. A synthetic 
photon spectra (solid line) generated from an electron power- 
law with £’, 0=20 ke'V and 6=4 (dashed curve in Fig. was 
fitted with a broken power-law (dashed line). The power-law 
left of the break energy e,o had a fixed spectral index of 1.5. 
For this example, fitting yields e,o=9 ke'V and a photon spectral 
index right of this break of y=2.91 (» 5-1). 

3. Known and unknown errors in computing 
non-thermai energies. 

So far, the emphasis has been on the difference due to the use 
of accurate cross-sections (both bremsstrahlung and Coulomb 
losses) or their NR approximations. In this section, other ap¬ 
proximations and factors are presented that also influence the 
photon spectrum and thus are possible sources of error for the 
derived electron energy distribution. Other effects on the pho¬ 
ton spectrum are non-uniform target ionization (Brown 1973) 
and the albedo effect (Compton back-scattering on the photo¬ 
sphere, Bai and Ramaty, 1978). They complicate the observed 
photon spectra by the simple fact that the amplitude of their 


effects varies with the energy. Kontar et al. (2002) for non- 
uniform target ionization and Alexander & Brown (2002) for 
the albedo effect have provided corrections. We have used their 
formulas to compute correction factors in the numerical exam¬ 
ples given below. 

Finally, the possible presence of some high-energy cutoff 
(or break) in the injected electron distribution affects the pho¬ 
ton spectrum at lower energies. The effect has been simulated 
numerically. A sharp electron high-energy cutoff at energy Eh, 
would lead to a noticeable deviation from power-law behavior 
starting already at photon energies above ~ where a pro¬ 
nounced rollover in the spectrum should occur. This was not 
observed in our selection of (mostly M-class) flares, at least 
below 35 keV, the upper limit of the fitting interval used in our 
data analysis. 

3.1. Some numerical examples 

How far from the truth are we if, from an observed photon 
power-law in the 10-35 ke'V band, we derive its electron power- 
law energy distribution using the Brown (1971) thick-target 
model (i.e. perfect power-law at all energies, non-relativistic 
cross-sections, cold target, no albedo, uniform 100 % ioniza¬ 
tion)? To estimate the effect of the different correction fac¬ 
tors, we have computed some numerical examples. Synthetic 
photon spectra were computed from ideal power-law electron 
distributions with index 6, using the relativistic cross-sections, 
and other effects (high-energy cutoffs, ionization, and albedo: 
See Table^. These synthetic photon spectra were fitted in the 
10-35 keV band (1-keV bins) with photon power-laws. Photon 
spectral indices, 7 , and fluxes at 50 ke'V, £ 50 , were determined 
from the fits. Using the Brown (1971) model, approximations 
to the original electron spectral indices (Eq.|3l and normaliza¬ 
tion factors (Eq.|2j can then be determined, from which non- 
thermal powers can be computed using Eqs. 0 or (do) (a 10 
keV cutoff or turnover energy was assumed here: This arbitrary 
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Type of correction 

Quantity 

6=3 

6=5 

6=7 

Haug, Bethe-Bloch 

Ay 

-t-0.1 

+0.1 

+0.1 


/so 

0.89 

0.91 

0.94 


fco 

1.07 

1.16 

1.23 


fro 

1.02 

1.14 

1.22 

Haug, Bethe-Bloch, 

Ay 

-0.2 

0.0 

+0.1 

Albedo 

/5O 

1.54 

1.36 

1.25 


fco 

0.96 

1.285 

1.40 


fro 

1.09 

1.29 

1.40 

Haug, Bethe-Bloch, 

Ay 

+2.1 

+ 1.2 

+0.7 

High-E cutoff 50 keV 

/so 

0.02 

0.2 

0.47 


fco 

2.64 

2.20 

1.71 


fro 

1.45 

1.95 

1.65 

Haug, Bethe-Bloch, 

Ay 

+0.35 

+0.1 

+0.1 

High-E cutoff 300 keV 

fso 

0.50 

0.90 

0.94 


fco 

1.23 

1.17 

1.23 


fro 

1.02 

1.15 

1.22 

Haug, Bethe-Bloch, 

Ay 

0.0 

0.90 

0.0 

E. = 5 keV 

/so 

2.54 

2.67 

2.84 


fco 

2.81 

2.74 

2.63 


fro 

2.74 

2.74 

2.64 

Haug, Bethe-Bloch, 

Ay 

-0.2 

-0.4 

-0.3 

E. = 25 keV 

/5O 

2.56 

2.30 

1.86 


fco 

1.51 

1.02 

1.01 


fro 

1.74 

1.09 

1.03 

Haug, Bethe-Bloch, 

Ay 

-0.2 

0.0 

+0.1 

£. = 100 keV 

/so 

1.74 

1.02 

0.96 


fco 

0.95 

1.10 

1.23 


fro 

1.12 

1.10 

1.22 

Haug, Bethe-Bloch, 

Ay 

-0.4 

-0.5 

-0.4 

E. = 25 keV, 

/5O 

4.42 

3.45 

2.46 

Albedo 

fco 

1.24 

1.12 

1.15 


fro 

1.84 

1.22 

1.18 


Table 1. Differences in photon spectral indices and normaliza¬ 
tion factor (flux at 50 keV), as well as computed non-thermal 
power (cutoff and turnover cases), between reality and assum¬ 
ing the Brown (1971) model. See text for details. in the 
context of the non-uniform target ionization model, is the min¬ 
imum initial energy an electron requires to reach the neutral 
parts of the chromosphere (same as in Brown 1973; Kontar et 
al. 2002 ). 


value does not greatly change the generality of the problem). 
We call these quantities dapprox, Ae^approx, PZwfp P‘lumover 
Ay (Table [0 is the difference between the spectral index y of 
the synthetic spectrum and d - 1 , the photon spectral index that 
would have been obtained with the Brown (1971) model (Eq. 
1^. Similarly, /50 is the ratio of the photon fluxes at 50 keV 
of the synthetic spectra with those derived theoretically using 
the Brown (1971) model (Eq.EJ. If P^t/f PturLer are the 
real non-thermal powers (derived from Eqs (|3 and (IIOH . the 
factors fco and fro'- 


fco - 


preal 

^cutoff 

j^approx 

^cutoff 


( 11 ) 


fro - 


preal 

^turnover 

j^approx 

^turnover 


( 12 ) 



Fig. 7. The photon turnover energy, eto (in keV), as obtained 
by htting synthetic spectra with a weak EM - 0.03 x 10"^® 
cm“^, T - \ keV thermal component, and a 6-5, A 50 = 10^^ 
electrons s ' keV ' non-thermal power-law distribution with 
different turnover energy Eto- 


are the corrections that must be applied to the rough 
papprox estimates. Table [I] indicates that the Ay - y - 

cutoj j,turnover ' • 

(d - 1 ) = Sapprox - 5 is usually slightly positive, i.e. applying the 
Brown (1971) model to observed photon power-law spectra to 
determine the original electron power-law spectral index gen¬ 
erally slightly underestimates that electron spectral index. In 
the 10-35 keV range, for usually observed spectral hardnesses 
{6 > 4), both corrections to the non-thermal power for rela¬ 
tivistic effects and albedo are of the same importance, about 
15-20%. The correction due to non-uniform target ionization is 
usually more important, particularly for low 6 and low E, {E, 
is the initial energy that electrons need in order to penetrate 
into the unionized chromosphere). Some of the effects, when 
combined, might cancel each other out, and all depend on the 
spectral index 6: General error estimates from each effect (or 
the sum of them) can only be done for a certain energy band of 
observation and if an approximate spectral index is known. The 
fro-co corrections on non-thermal power exceed 50% only in 
extreme cases. The biggest uncertainty comes by far from the 
low-energy cutoff or turnover energy: the non-thermal power 
going as The addition of albedo and/or non-uniform tar¬ 

get ionization effects combined with the fact that spectral htting 
is somewhat model-dependant may displace this Eco,to- 

3.2. Finding the iow-energy cutoff or turnover 

This paragraph assumes the turnover model, but conclusions 
are qualitatively the same for the cutoff model. As shown in 
Figs-HIH] and 1 ^ the presence of a thermal component in the 
photon spectrum complicates the determination of Eto, the elec¬ 
tron turnover energy. An unambiguous relationship between 
and Eto cannot always be established. Making several httings 
with different T, EM, 6, A 50 , and Eto yielded the following rule 
of thumb: e,o, and the Eto derived from it, seem reliable only 
when eto > ^thnnth . where ethnnth is the energy where the ther- 
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Fig. 8. As in Fig.0 T - \ keV but with a very strong emission 
{EM - 10'*® cm“^) for the thermal component. For Eto <18 
keV, the eto fitting parameter may take any value (below 15 
keV), and still yield a very good fit. 



Fig. 9. Synthetic spectrum obtained with Eto - 20 keV, other 
parameters as in Fig. The Haug and Bethe-Bloch cross- 
sections have been used. 

mal and non-thermal components of the spectrum intersect. In 
a case such as depicted by Figs. Elandli if the fitted is below 
15 keV, then only a upper value of ~18 keV may be assigned 
for Eto, leading to a lower boundary for the non-thermal power. 

4. RHESSI flare observations, data analysis 

4.1. Flare selection criteria 

In the remaining part of this paper, the effects of the corrections 
and improvements on the computation of electron collisions 
and bremsstrahlung emission are studied on real data. The di¬ 
rection is now reversed: from the photon spectrum the electron 
energy distribution is derived and the total non-thermal energy 
is estimated. The results of the previous section serve to esti¬ 
mate errors. 

Flares have been selected from the first 18 months of 
RHESSI observations, using the following selection criteria: 


Flare 

Time intervals 
studied (approx.) 

pos. 

[”] 

2002/04/09 12:59 

12:57:15-13:02:20 

-569,405 

2002/06/01 03:53 A 
2002/06/01 03:53 ABC 

03:53:40-03:54:35 

03:53:40-03:58:50 

-423,-303 

2002/07/11 14:18 A 
2002/07/11 14:18 AB 

14:17:15-14:18:10 

14:17:15-14:19:50 

-791,281 

2002/08/22 01:52 A 
2002/08/22 01:52 B 

01:49:25-01:50:10 

01:50:10-01:55:20 

816,-272 

2002/09/08 01:39 

01:37:10-01:40:20 

-911,-205 

2002/10/05 22:50 

22:50:00-22:52:45 

-558,72 

2002/11/10 03:11 

03:07:00-03:15:00 

592,-240 

2002/11/14 11:09 A 
2002/11/14 11:09 A 

11:09:15-11:10:35 

11:09:25-11:09:57 

-887,-262 

2003/06/10 02:52 ABC 
2003/06/10 02:52 B 
2003/06/10 02:52 C 

02:51:15-02:53:50 

02:51:45-02:52:20 

02:52:20-02:53:25 

561,185 


Table 2. Flares, the time intervals that were used, and their an¬ 
gular offsets from Sun center. 


- For simplicity’s sake from a data analysis point of view, 
flares (or portions of flares) with the same attenuator states 
throughout (including background time), no decimation 
and no pile-up were taken. 

- Flares had to have two foot points, in order to determine 
a loop volume. HEDC' images (7” resolution, in different 
energy bands) were used to determine this. 

- Significant HXR flux above 25 keV was required. 

- Only flares above C5.0 GOES X-ray level were selected. 

From the many cases 9 flares have been selected. They were 
relatively simple, but some have more than one HXR >25 keV 
episode, in which case the peaks were labeled chronologically 
A, B, C,... All selected flares turned out to have an attenuator 
state of 1 (thin shutter in). Table|2]lists them. 

4.2. Extracting the thermal and non-thermal flare 
energies 

4.2.1. Spectral fitting 

Using the SPEX software package of the RHESSI standard 
analysis tools, spectral models composed of a thermal com¬ 
ponent and a broken power-law are fitted to RHESSI spectra. 
The low-energy power-law has a fixed spectral index of 1.5. 
This somewhat arbitrary value is an average approximation of 
photon spectral indices at photon energies below the turnover 
energy. The fitting was done for time intervals varying between 
2 to 5 RHESSI spin periods (~4 s), in the 6-35 keV band. 
This band was chosen because lower energies could depend too 
much on the model used and the accuracy of the instrument’s 
spectral response matrix. At higher energies, a spectral break 
may be present (e.g. Fig. [Hi. a high-energy spectral break is 
most likely due to a break in the original injected electron dis¬ 
tribution (Miller 1998). For simplicity, the details of the spec- 

' http://www.hedc.ethz.ch 
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Fig. 10. RHESSI photon spectrum, taken near the beginning of 
HXR emission of the 2002/09/08 Ml .6 flare. 


trum above this spectral break were omitted as their influence 
on the non-thermal power is negligible. 

The five parameters retrieved are the temperature T and 
emission measure EM of the thermal component, the spectral 
index y, normalization factor Fso, and the turnover photon en¬ 
ergy e,o of the power-law component. Only time intervals with 
significant HXR flux above 25 keV were studied. This require¬ 
ment allows better accuracy in determination of non-thermal 
energies. It also has the effect of dividing the flare into episodes 
of large energy input for time intervals that are short, and hence 
with smaller thermal energy losses. 

During fitting, high values for ew were started with, tending 
to yield an upper boundary for this value in cases where it was 
not well-defined, as happens when the photon turnover occurs 
near or within the thermal part of the spectrum (which leads to 
a lower limit for the computed non-thermal energy). 

4.2.2. Computing non-thermal energies 

In the following, the turnover model is considered. Every time 
interval yields the following fitting parameters: y, Fso, Sro, T, 
and EM. To compute non-thermal power or energy from Eq. 
(GOll, the following must be determined, d, Ag, and Efo. 

Ag and d: To determine the correct Ag and 6 to use, the follow¬ 
ing method has been used, for every data point: 


- Erom y' and A', and again using Eqs. (|2i and 0, an in¬ 
jected electron spectral index d' and normalization factor 
A' are determined. 

- We define: 


Ad = 

^ ^approx 

(13) 


K 

(14) 

Pe — 



If one uses the Brown (1971) method (i.e. Eqs|2and0 to 
determine the injected electron spectral characteristics from 
an observed photon spectrum, then the electron spectral in¬ 
dex will have been overestimated by A(5 (underestimated if 
Ad < 0 ), and the normalization constant by pg (see also 
TableHJ. 

- Ad and pg depend mostly on the photon spectral index, and 
vary slowly with it. Assuming that Ad is small and pg near 
unity, so that the relationship with y is linear, the increase 
in spectral index from d (the real original injection spec¬ 
tral index we are looking for) to dappmx is about the same 
as from dappmx to d', i.e. Ad. Similarly, the increase in 
the normalization constant is similar from Ag to Ag^appmx as 
from Ag^approx to A', i.e. pg. Hence, the d and Ag to be used 
in Eq. (HI can be approximated by: 


d ^ dappmx - Ad 

(15) 

. ^e,approx 

Ag ~ 

(16) 


Pe 


Non-uniform target ionization effects have been neglected 
here: This correction seems unnecessary in light of the fact 
that no energy break in the relevant energy band (< 35 keV) 
for fitting were observed. This could be due to the fact that 
expected features from non-uniform target ionization in the 
observed photon spectrum might go unnoticed (masking by 
thermal emission at the low energies, by count statistics at 
the high energies): In this case, the non-thermal energies 
could be overestimated by up to a factor ~2.8 (Table 0 
for small £* (below ~5 keV). Assuming £, >25 keV, and 
electron spectral index d > 3.5, the error is at most +20%, 
and progressively less as £, and/or d increase. 

Eto'. To And the correct Eto, the following has been done: 

- Using, d and Ag as determined above, as well as the T and 
EM fitting parameters, a graph such as those presented in 
Figs.0and|3is generated. 

- As explained previously, one can find E,o from e^, or at 
least an upper limit for it. 


Erom the y and Fso fitting parameters, Eqs. 0 and 0 
(with Ag - 50^ ■ Fso), are used to determine an approximate 
injected electron spectral index dappmx and normalization 

constant Ag^appmx- 

Using these dappmx and Agappmx, a synthetic photon spec¬ 
trum using relativistic cross-sections and the albedo cor¬ 
rection is computed. The resulting photon spectrum is now 
fitted by a power-law the same way as the observed data. 
This yields a spectral index y' and normalization constant 
A' 


With d, Ag, and E,o, the non-thermal power and energy can 
be computed. 

4.2.3. Computing thermal energies 

Thermal energies are computed using: 

E,h = 'iksT xjEM -V-f (17) 

Equal electron and ion temperatures and a unity filling factor 
/ were assumed throughout. The assumption of a near-unity 
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Flare 

Class 

Flare volume 
[xlO-’ cra^] 
Min Max 

2002/04/09 12:59 

M1.2 

1.0 

4.7 

2002/06/01 03:53 

M1.6 

0.48 

3.6 

2002/07/11 14:18 

Ml.O 

1.3 

4.2 

2002/08/22 01:52 

M7.8 

1.9 

9.6 

2002/09/08 01:39 

M1.6 

0.55 

9.7 

2002/10/05 22:50 

C6.8 

0.12 

4.2 

2002/11/10 03:11 

M2.6 

0.66 

9.3 

2002/11/14 11:09 

C5.9 

0.38 

1.7 

2003/06/10 02:52 

M2.2 

2.5 

3.8 


Table 3. Flare volumes, as determined (Min) by RHESSI 6-8 
keV imagery, or (Max) from non-thermal energies (see text). 


filling factor is supported by e.g. Dere (1982) or Takahashi 
& Watanabe (2000). On the other hand, there seem to be ev¬ 
idence that filling factors as low as 10“^ or even 10“^ exist 
(Cargill & Klimchuk 1997), and ~0.1 is an oft reported value. 
Two methods were used to determine the flare volume V, both 
using RHESSI images done with the CLEAN algorithm, and 
with different sets of collimators (1-7, 2-7 and 3-7). The first 
one consisted in making an image in the 6-8 keV band (all 
flares used in this study had a visible Fe-Ni line complex above 
the free-free continuum), estimating the flare area A, and us¬ 
ing V - as the flare volume. The second one consisted in 
looking at images made at non-thermal energies (usually 25-50 
keV), estimating the size of foot points by fitting 2-D elliptical 
gaussians, then deconvolving for the CLEAN beam. A volume 
is computed from the sizes and the distance between the foot 
points, assuming a perfect arc-shaped loop. The first method 
yielded a lower value for V and the thermal energy content, 
while the second one provided an upper value for both. The 
range of values for V easily reached an order of magnitude. 
Considering RHESSFs dynamic range, a loop 10 times larger 
than a smaller one may be invisible: Assuming both loops have 
the same plasma content, the surface brightness of the larger 
one is about ~50 times lower than a that of the smaller one 
(cf. AnnendixlAl. This justifies considering the second method. 
Tables 13 and 0] list the flare volumes and the thermal energies 
derived from them. The thermal energy increases between the 
start and the end of HXR flux >25keV were considered: 

AEih - Eth^HXRend - Efh^HXRstart ( 18 ) 

Eigure[nis an example of the acquired data and some de¬ 
rived quantities for the 2002/11/10 flare. Other flares studied 
yielded qualitatively similar behaviours, and are not displayed 
here. 

5. Results and discussion 

As displayed in Eig. in the peak in the 50 keV flux almost 
always coincides with a dip in the spectral index of the pho¬ 
ton power-law. This common behaviour is thought to be a con¬ 
sequence of the acceleration process (Grigis & Benz, 2004). 
As observed in all our flares, the emission measure EM (both 
RHESSI- and GOES-derived) increases during the flare. The 




03:08 03:10 03;12 03:14 


Fig. 11. Data for the 2002/11/10 flare, during the main HXR 
peak. From top to bottom: (1) Light curves in RHESSI counts 
at low, intermediate, and high photon energies; (2) emission 
measure of thermal plasma; (3) temperature of thermal plasma; 
(4) calibrated photon flux at 50 keV; (5) turnover electron en¬ 
ergy (upper limit); (6) photon spectral index; (7) non-thermal 
power; (8) cumulative non-thermal energy and minimal ther¬ 
mal energy as derived from RHESSI or GOES temperature and 
emission measures. 

RHESSTderived temperatures are usually above the GOES- 
derived ones, while the RHESSTderived emission measures 
are below the GOES-derived ones. However, the T sjEM prod¬ 
uct is similar most of the time (as can be deduced from the 
thermal energies of Tabled}- The RHESSTderived temperature 
initially decreases rapidly, then stabilizes during the rest of the 
time when significant HXR flux > 25 keV is present. This ini¬ 
tial decrease in RHESSTderived temperature is not always ob¬ 
served, and is never present in the GOES-derived temperatures. 
The RHESSI spectral fittings were often less reliable at those 
early times, leading the authors to believe that the RHESST 
derived T and EM values are not reliable at those early times. 
At the later times, thermal energies derived from both RHESSI 
and GOES T and EM yield similar values. It might be argued 
that the isothermal bremsstrahlung spectrum might not always 
be the best model for fitting the thermal component of X-ray 
spectra, and that fitting a multi-thermal model, or a full dif¬ 
ferential emission measure distribution would be more proper, 
although practically more difficult. 
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Flare 

Class 

dt[s] 

EnthJO 
xlO^® ergs 

^^th,RHESSI,mini max 

xlO^” 

^EtfJi^QQES jnin ! max 

ergs 

2002/04/09 12:59 

M1.2 

309.4 

2.0 

0.49/1.06 

0.48/1.04 

2002/06/01 03:53 A 

M1.6 

98.2 

0.40 

0.19/0.52 

0.25/0.68 

2002/06/01 03:53 ABC 

M1.6 

307.0 

0.47 

0.58/1.59 

0.44/1.21 

2002/07/11 14:18 A 

Ml.O 

48.1 

0.05 

0.41/0.74 

0.37/0.66 

2002/07/11 14:18 AB 

Ml.O 

120.3 

0.48 

0.48/0.86 

0.44/0.79 

2002/08/22 01:52 A 

M7.8 

37.2 

0.09 

0.04/0.09 

0.048/0.11 

2002/08/22 01:52 B 

M7.8 

298.6 

21.6 

1.74/3.9 

1.53/3.4 

2002/09/08 01:39 

M1.6 

162.0 

6.9 

0.54/2.3 

0.46/1.9 

2002/10/05 22:50 

C6.8 

158.7 

0.12 

0.16/0.95 

0.21/1.24 

2002/11/10 03:11 

M2.6 

468.1 

4.5 

0.58/2.2 

0.55/2.1 

2002/11/14 11:09 A 

C5.9 

72.4 

0.20 

0.18/0.38 

0.21/0.44 

2002/11/14 11:09 A 

C5.9 

32.2 

2.1 

1.0/2.1 

1.1/2.3 

2003/06/10 02:52 ABC 

M2.2 

125.9 

6.0 

0.62/0.76 

0.48/0.59 

2003/06/10 02:52 B 

M2.2 

33.6 

0.39 

0.19/0.23 

0.15/0.18 


Table 4. Time intervals of significant HXR flux >25 keV, total corrected non-thermal energy for either the cutojf or the turnover 
models, ratio of non-thermal energy over the thermal energy increase during that time interval, with the thermal energy increase 
computed using either RHESSI or GOES T, EM data. The numbers in parenthesis are the range of values, due to the uncertainty 
in the thermal volume. 


The dashed line in the fifth plot of Fig. IIII corresponds to 
times where the photon turnover energy 6,0 is clearly above 
the thermal part, hence yielding a reliable value for the corre¬ 
sponding electron turnover energy Eto- Later, the thermal part 
becomes so important that only an upper value for Eto may be 
determined. The turnover energy Eto does not seem to change 
substantially during the main HXR phase, and increases to a 
higher value later in the flare (similar to the 2002 July 23, 
Holman et al. 2003). The derived time of peak non-thermal 
power does not exactly coincide with the time of peak pho¬ 
ton emission at 50 keV for this flare. This may not be real, as 
the electron turnover energy Eto is only an upper limit at those 
times. 

The turnover model yields non-thermal energies typically 
only 10 to 30 % higher than the cutoff model. This stems from 
the fact that fitting our double photon power-law on the photon 
spectra from a turnover electron model always yields a larger 
photon turnover eto than with a cutoff model. This translated to 
a higher Eto than Eco, leading to turnover-model non-thermal 
power only slightly higher than the cutoff-model non-thermal 
power. 

The total non-thermal energy and thermal energy increase 
for the studied flares are summarized in Tabled] 

Table 13 lists the non-thermal to thermal energy ratios, for 
different durations of the HXR peak. The non-thermal energies 
are lower limits, the non-thermal to thermal ratios are hence 
also lower limits. Ratios obtained using the minimal flare vol¬ 
ume are arguably closest to the truth, most notably because of 
filling factor considerations. It can be noted that short HXR 
peaks lead to ratios of ~1.5, whereas longer-duration peaks 
lead to higher ratios: ~6. This is expected, as radiative cool¬ 
ing or heat conduction (the second being most likely the domi¬ 
nant loss mechanism: Cargill 1994, Porter and Klimchuk 1995) 
tend to lower the thermal energy content, thereby increasing the 
ratio. Taking time intervals ending well after the main HXR 



0.01 0,10 1.00 10,00 100.00 
E„„ [xio” ergs] 



0.01 0.10 1.00 10.00 100.00 
Enih [x1 0^° ergs] 


Fig. 12. Log-log plot of non-thermal vs. thermal energies for 
all flares in Table 0] The solid line is a linear fitting, yielding 
constant a and slope b. The dashed line is also a linear fitting, 
using the bisector method. 


peaks, such as the SXR peak, would tend to lower the thermal 
energy, and to increase the non-thermal/thermal ratio. 

All data points have also been plotted on Fig. El Linear 
fitting with the bisector method (Isobe et al. 1990), more rele¬ 
vant in cases where variables are truly independent, yields the 
following relation: 

^Eth ~ ( 19 ) 

This empirical relation may simply state that the thermal 
energy increase AEth does not increase as fast as the cumulative 
non-thermal energy, due to losses. 
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Flare 

volume 

T,EM 

source 

Ratio of non-thermal to thermal 
All HXR peaks t<75s 75s<t<200s 

energies 

t>75s 

t>200s 

min. 

RHESSI 

4.2±4.5 

1.5±0.9 

5.3±5.6 

5.7±5.0 

6.3±5.0 


GOES 

4.7±5.4 

1.4±1.0 

6.1±7.0 

6.5±6.0 

6.9±5.6 

max. 

RHESSI 

1.6±1.6 

0.7±0.4 

1.7±1.8 

2.1±1.9 

2.4±2.2 


GOES 

1.7±2.0 

0.7±0.4 

2.1±2.3 

2.3±2.3 

2.7±2.5 


Table 5. Ratio of non-thermal to thermal energies, for different duration t of HXR emission. 


6. Conclusions 

The investigation of the various cross-sections approximations 
(relativistic vs. non-relativistic) has shown that they played a 
relatively minor role for our energy estimation, under usual 
flare conditions. They are important, however, for spectral in¬ 
versions and accurate derivation of the injected electron spec¬ 
trum. The determination of the low-energy cutoff or turnover, 
the largest source of error for non-thermal energies, seems re¬ 
liable in some flares at the beginning of HXR emission. With 
the procedure that was used in this paper, determining the non- 
thermal power assuming a turnover model yielded usually only 
slightly more (~20-25% on the average) energy than if assum¬ 
ing a cutoff model. 

The potentially largest sources of error for thermal energy 
estimations are the uncertainty in Ailing factors, and the fact 
that energy losses play a substantial role over long accumula¬ 
tion times. 

The non-thermal energies and the thermal energies in a flare 
seem to be of the same order of magnitude, at least during the 
initial stages. Later on, radiative cooling or heat conduction 
lower the thermal energy content, thereby increasing the ratio. 
The ratio of non-thermal over thermal energies given in Table 
13 are lower limits: Both lower low-energy cutoffs/turnovers 
or Ailing factors might increase it substantially. A bigger than 
unity ratio was expected because evaporation includes heating 
and expansion (Benz & Krucker 1999), and the non-thermal 
to thermal energy conversion might not be 100% effective 
(some of the energy might rapidly be lost radiatively by cooler 
plasma). 
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Appendix A: Surface brightness of two volumes 
with the same plasma contents 

The surface brightness f of a feature is proportional to EM I A, 
where EM is the emission measure (EM - rP'V), and A the 
feature’s area. Let A = Assuming a constant plasma con¬ 
tent {nV) and temperature, F is hence proportional to 
Two features possessing the same plasma content, but of dif¬ 
fering sizes, may have a widely differing surface brightness: if 
the first feature is 10 times smaller in volume, it will have 46.4 
times larger surface brightness! The total thermal X-ray flux 
from the smaller feature will also be 10 times larger. 







